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Abstract. A self-consistent computational scheme is presented for one dimensional 
(ID) and two dimensional (2D) metamaterial systems with gain incorporated into 
the nanostructures. The gain is described by a generic four-level system. The loss 
compensation and the lasing behavior of the metamaterial system with gain are studied. 
A critical pumping rate exists for compensating the losses of the metamaterial. There 
exists a wide range of input signals where the composite system behaves linearly. 
Nonlinearities arise for stronger signals due to gain depiction. The retrieved effective 
parameters are presented for one layer of gain embedded in two layers of Lorentz 
dielectric rods and split ring resonators with two different gain inclusions: (1) gain is 
embedded in the gaps only and (2) gain is surrounding the SRR. When the pumping 
rate increases, there is a critical pumping rate that the metamaterial system starts 
lasing. 
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1. Introduction 

The field of metamaterials [U [2] is driven by fascinating and far-reaching theoretical 
visions such as, e.g., perfect lenses [3], invisibility cloaking [HE], and enhanced optical 
nonlinearities [6]. This emerging field has seen spectacular experimental progress in 
recent years [U [2]. Yet, losses are orders of magnitude too large for the envisioned 
applications. Achieving such reduction by further design optimization appears to be 
out of reach. Thus, incorporation of active media (gain) might reduce the losses. The 
procedure would be to simply inject an electrical current into the active medium, leading 
to gain and hence to compensation of the losses. However, experiments on such intricate 
active nanostructures do need guidance by theory via self-consistent calculations (using 
the semi-classical theory of lasing) for realistic gain materials that can be incorporated 
into or close to dispersive media to reduce the losses at THz or optical frequencies. The 
need for self- consistent calculations stems from the fact that increasing the gain in the 
metamaterial, the metamaterial properties change, in turn changes the coupling to the 
gain medium until a steady-state is reached. A specific geometry to overcome the severe 
loss problem of optical metamaterials and to enable bulk metamaterials with negative 
magnetic and electric response and controllable dispersion at optical frequencies is to 
interleave active, optically pumped gain material layers with the passive metamaterial 
lattice. 

For reference, the best fabricated negative-index material operating at around 
1.4 /iin wavelength [7] has shown a figure of merit, FOM = — Re(n)/Im(n) « 3, where 
n is the effective refractive index. This experimental result is equivalent to an absolute 
absorption coefficient of a = 3 x 10 4 cm _1 , which is even larger than the absorption of 
typical direct-gap semiconductors such as, e.g., GaAs (where a = 10 4 cm _1 ). So it looks 
difficult to compensate the losses with this simple type of analysis, which assumes that 
the bulk gain coefficient is needed. However, the effective gain coefficient, derived from 
self-consistent microscopic calculations, is a more appropriate measure of the combined 
system of metamaterial and gain. Due to pronounced local-field enhancement effects 
in the spatial vicinity of the dispersive metamaterial, the effective gain coefficient can 
be substantially larger than its bulk counterpart. While early models [H [9j [10l [HJ [12] 
using simplified gain-mechanisms such as explicitly forcing negative imaginary parts of 
the local gain material's response function produce unrealistic strictly linear gain, our 
self-consistent approach presented below allows for determining the range of parameters 
for which one can realistically expect linear amplification and linear loss compensation 
in the metamaterial [13] . To fully understand the coupled metamaterial-gain system, we 
have to deal with time- dependent wave equations in metamaterial systems by coupling 
Maxwell's equations with the rate equations of electron populations describing a multi- 
level gain system in semi-classical theory [14J. 

This paper aims to apply a detailed computational model to the problem of 
metamaterials with gain. In section [2], we present the semi-classical theory of lasing 
and describe in detail the computational approach. In section [31 we verify that our 
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Figure 1. (Color online) Schematic of the four-level atomic system model. 

code agrees well with simple soluble models (gain material only). In addition, our code 
is applied to ID superlattice of gain and negative index layers. Next, a 2D problem 
is considered, which is a square lattice of Lorentz dielectric cylinders with layers of 
gain material. Finally, a 2D split ring resonator (SRR) with gain material inclusions 
is considered. Gain can compensate the losses and lasing (spasing) is achieved in our 
numerical simulations. In section HJ we present our conclusions. 



2. Theoretical and numerical model 

The gain atoms are embedded in host medium and described by a generic four-level 
system, as shown in figure HJ All quantities including the fields and occupation numbers 
are tracked at each point in space and take into account energy exchange between 
gain atoms and fields, external pumping and non-radiative decays [14J. Electrons are 
pumped by an external mechanism from the ground state level (Ao) to the third level 
(A3). After a short lifetime T32, they quickly relax into the metastable second level 
(iV 2 ). The second level (A 2 ) and the first level (Aq) are called as the upper and lower 
lasing levels, respectively. Electrons can transfer both radiatively (spontaneous and 
stimulated emissions) and non-radiatively from the upper to the lower lasing level. At 
last, they transfer quickly and non-radiatively from the first level (Aq) to the ground 
state level (A ). The energies of ground state and the third level are E and E 3 . In 
optical pumping mechanism, electrons are raised from the ground state level (A ) to 
the third level (A3) by an external electromagnetic wave with the pumping frequency 
ojb = (E 3 — Eq)/H, which is chosen to be 4tt x 10 14 Hz in our simulations. The 
local intensity of the pumping EM wave varies with the position and determines the 
pumping rate at each point. The lifetimes and energies of the upper and lower lasing 
levels are r 2 i, E 2 and r 10 , E 1 , respectively. The center frequency of the radiation is 
u a = {E 2 — Ei)/h, which is chosen to be 2it x 10 14 Hz. The parameters T32, T21 and 
t 10 are chosen 5 x 10~ 14 s, 5 x 10~ 12 s and 5 x 10~ 14 s, respectively. The total electron 
density, N (t = 0) = N (t) + Ni(t) + N 2 {t) + A 3 (t) = 5.0 x 10 23 /m 3 . 
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The time-dependent Maxwell equations are given by 

V x E = -dB/dt 

1 

V x H = ee dE/dt + dP/dt, 

where B = yUyU H and P = Y2i= a b * s ^ ne electric polarization density of the gain 
material. (P a is the induced electric polarization density on the atomic transition 
between the upper (N 2 ) and lower (N{) lasing levels, and P b is between the ground 
state level (No) and the third level (N 3 ).) The induced electric polarizations behave as 
harmonic oscillators and couple to the local E field, which is propagated by Maxwell 
equations. The polarization density Pj(r, t) obeys locally the following equation of 
motion [H] 

d 2 P(t) dP(t) 

^f 1 + T.^P + w?P 4 (t) = -a t AN(t)E(t) (i = a,b), (2) 

where Tj is the linewidth of the atomic transition u>i, &i is the coupling strength of Pj to 
the electric field, and AiV a (r, t) = N 2 (r, t) - N ± (r, t) and AiV b (r, t) = JV 3 (r, t) - iVo(r, t) 
are the population inversions that drive the polarizations. In our simulations, T a is 
chosen to be equal to 2ir x 5 x 10 12 Hz or 2ir x 20 x 10 12 Hz and T b is equal to 
2n x 10 x 10 12 Hz. The values for a a and a b are taken to be 10~ 4 C 2 /kg and 5 x 10~ 6 C 2 /kg, 
respectively. From ([2]), it can be easily derived [14J that the atomic response of gain 
atoms has a Lorentzian lineshape and is homogeneously broadened. The occupation 
numbers at each spatial point vary according to the following rate equations, 
dN 3 
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where -^-E • ^ (i = a, 6) is the induced radiation rate or excitation rate depending on 
its sign. 

Instead of using an external EM wave to optically pump electrons from the ground 
state level (N ) to the third level (iV" 3 ), we can simplify this process in fl3aj) and (j3dj) by 
pumping electrons with a homogeneous pumping rate r pump , which is proportional to 
the optical pumping intensity in an experiment. This simplification is valid only if the 
gain slab is thin and the gain of the laser is low, because the real pumping rate depends 
on the local optical intensity and should be a function of position. We'll discuss this in 
more detail in section 13.11 Based on this simplification, we can have the rate equations 
as follows, 

dN 3 N 3 

= r pump iVo , (4a) 

ot r 32 
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Figure 2. (Color online) Schematic of gain material slab (shown in orange). The slab 
width w takes different values in the cases we have examined. 



9N 2 N 3 1 ^ dP a N 2 
dt r 32 Ak^a at r 2 i 

^ = ^-^E-^-^, (4c) 
<9t r 2 i /k^ a <9i no' 

^ = ^ - r pump iv , (4d) 

at r 10 

Correspondingly, we only need to consider the electric polarization density P (r, t) on 
the atomic transition between N 2 and N\ in ([1]) and (J2]). 

In order to solve the behavior of the gain materials in the electromagnetic fields 
numerically, the finite-difference time-domain (FDTD) method is utilized [151 HSl [T7] . 
At the left and right ends of the computational space, perfect matched layers (PML) are 
used to impose the absorbing boundary condition (ABC). In the FDTD calculations, 
both the space and time are discretized into small steps compared to the characteristic 
space and time periods. In our simulations presented below, the discrete time and space 
steps are chosen to be At = 1.67 x 10~ 17 s and Ax = 1.0 x 10~ 8 m in section 13.11 
and section 13.21 At = 8.33 x lCT 18 s and Ax = 5.0 x 10 _9 m in section 13. 3[ and 
At = 8.33 x 10~ 19 s and Ax = 1.0 x 10~ 9 m in section I3T41 The initial condition is that 
all electrons are in the ground state and all electric, magnetic and polarization fields are 
zero. Then the electrons are pumped from N to iV 3 optically or with a homogeneous 
pumping rate r pump . The system begins to evolve according to the equations above. 



3. Examples for metamaterials incorporated with gain 

3.1. Gain material only 

To understand the lasing behavior of gain material, we first study a gain material slab 
surrounded by vacuum (shown in figure [2]). We generate a continuous wave (CW) at 
the frequency u b (200 THz) and let it propagate through the gain slab, and then we 
calculate the reflected and transmitted waves and implement the Fourier transforms to 
see if there is lasing and how much power is emitted around 100 THz — the emission 
frequency oj a between N\ and N 2 . First we start with a very low input power P in for 
the incident CW wave, but no lasing happens, then we increase the input power till it 
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Figure 3. (Color online) The transmitted waves and their corresponding Fourier 
transforms for different input powers, (a), (b) and (c) are the transmitted waves 
for input power P m = 79.6, 90.7 and 120.6 W/mm 2 , respectively, (d), (e) and (f) are 
same as (a), (b) and (c), respectively, but for the Fourier transforms of the transmitted 
waves. The gain slab width w = 100 nm and the bandwidth r a of the atomic transition 
between Ni and N2 is 5THz. 



reaches the lasing threshold, for which the system starts to have lasing and we can see 
a small peak at the emission frequency 100 THz in the Fourier transforms of reflected 
and transmitted waves, i.e., low power emitted around the emission frequency u a . If 
we keep increasing the input power, the peak will get higher and the emitted power 
will get larger. Figure [3] shows the transmitted waves and their corresponding Fourier 
transforms for three different input powers at the gain slab w = 100 nm. We can see 
there is no lasing (figure 3(a)) when the input power is low (P; n = 79.6 W/mm 2 ) and 
there is only one peak for the pumping frequency in its Fourier transform (figure 3(d) ). 
When the input power P in = 90.7 W/mm 2 , the system starts lasing (figure 3(b)) and 
a small peak appears at the frequency 100 THz (figure 3(e)). If we increase the input 
power to a higher value P in = 120.6 W/mm 2 , the lasing gets stronger (figure 3(c)) and 
the peak for the emission frequency gets higher (figure 3(f)), i.e., more power emitted 
around the emission frequency uj a . We have calculated the emitted power at the emission 
frequency u a versus the input power at the pumping frequency u>b for the same gain slab 
system. As shown in figure HI we can see that there is a sharp rise in the emission 
around P in ~ 90.7 W/mm 2 , which corresponds to the lasing threshold for this system. 
Below the threshold, there is no lasing. 

We also notice that the lasing time (the time when the system starts lasing) varies 
according to the input power. Figure [5] shows the detailed results for the lasing time 

100 nm, 250 nm and 500 nm. We can 



versus the input power with the slab width w 
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Figure 4. (Color online) The powers emitted at the emission frequency lu = u a 
(100 THz) for different input powers at the pumping frequency u> — uib (200 THz). All 
parameters of this system are same as figure [31 
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Figure 5. (Color online) The lasing times for different input powers at the pumping 
frequency lo^ (200 THz). The gain slab width w — 100 nm, 250 nm and 500 nm, 
respectively. All other parameters are same as figure |3j 

see the lasing time decreases as the input power increases because the system pumps the 
electrons at a higher rate from the ground state level to higher levels and then reaches 
the population inversion between Ni and N 2 in a shorter time. For a fixed input power, 
one can see the lasing time decreases as the gain slab width gets larger. This occurs 
because more input energy is absorbed and then converted to lasing by the wider gain 
slab. 

As the input wave propagates inside the gain slab, it will decay due to the absorption 
from the gain material at the pumping frequency u>b (see figure [6]). Thus the pumping 
rate, which is determined by the local input optical intensity, is inhomogeneous inside 
the gain slab. But for a thin gain material layer, the electric field of the input wave can 
be approximately treated as homogeneous, thus we can simplify the pumping process 
between No and N 3 by using a homogeneous pumping rate r pump . For a very thin gain 
slab w = 100 nm, simulations are done with a homogeneous pumping rate and the results 
for the power emitted around the emission frequency uo a versus the pumping rate are 
plotted in figure [71 Comparing with figure HJ where the electrons are optically pumped, 
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Figure 6. (Color online) The amplitude of the input EM wave inside the gain slab 
as a function of the position. The gain slab width w = 1000 nm and the input power 
Pin = 92.3 W/mm 2 . All other parameters are same as figure [3l 
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Figure 7. (Color online) The powers emitted at the emission frequency lu = u a 
(100 THz) for different pumping rates. The gain slab width w = iOOnm and the 
bandwidth of the atomic transition between N\ and iVjj is 5 THz. 



we can see they are very similar. For a fixed output power such as P out = 7.56 W/mm 2 , 
we can find the corresponding input power P in = 120.6 W/mm 2 in figure H] and the 
corresponding pumping rate r pump = 9.3 x 10 9 s _1 in figured Then we do simulations for 
optical pumping case with the input power P in = 120.6 W/mm 2 and for the homogeneous 
pumping rate case with the pumping rate r pump = 9.3 x 10 9 s _1 . The graphs of the 
occupation numbers as a function time are plotted in figure M for both cases. One can 
see that they are almost the same. This verifies that the homogeneous pumping rate 
simplification is valid for a thin gain slab. In our following simulations, we'll use this 
simplification because the gain slab widths in our structures are very thin (w <= 50 nm). 



3.2. Negative index material (NIM) embedded in layers of gain 

As the first simple model system, we will discuss a one- dimensional metamaterial system 
which consists of layers of negative index material (NIM) and gain material, as shown 
in figure [91 to see if we can compensate the losses of the metamaterials associated with 
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Figure 8. (Color online) The normalized occupation numbers as a function of time. 
The gain slab width w = 100 nm and the gain bandwidth of the atomic transition 
between ATj and N 2 is 5THz. (a) the electrons are optically pumped by an input 
EM wave with input power Pj n = 120.6 W/mm 2 and (b) the electrons are pumped 
with a homogeneous pumping rate T 



pump 



= 9.3 x 10 J s . Occupation numbers -/V , 
Ni, N 2 and 7V 3 are normalized by the total electron density Ni [Ni = N (t = 0) = 
5.0 x 10 23 /m 3 ]. 
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Figure 9. (Color online) The negative index material (blue) embedded in layers of 
gain material (orange). The number of layers, the permittivity and permeability of 
NIM are taken different values for different cases we have examined. The width for 
both NIM and gain material is w — 50 nm. The gain bandwidth is 5THz. 



the NIMs by the amplification provided by the gain material layers and how the system 
starts lasing. We first let a narrow band Gaussian pulse of a given amplitude go through 
the metamaterial without gain, and we calculate the transmitted signal emerging from 
the metamaterial system, which also has a Gaussian profile but its amplitude is much 
smaller than that of the incident pulse due to the losses of NIM layers. Then we introduce 
the gain into the system and start increasing the pumping rate. The amplitude of the 
transmitted signal gets larger and we can find a critical pumping rate, for which the 
transmitted pulse is of the same amplitude as the incident one. Since the gain material 
itself is nonlinear, we increase the amplitude of the incident Gaussian pulse for a fixed 
pumping rate to see how high we can go in the strength of the incident electric field 
and still have the full compensation of the losses, i.e., the transmitted signal equals the 
incident signal, independent on the signal strength. In this region we have compensated 
loss and still have linear response of the metamaterial. The shape of the transmitted 
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Figure 10. (Color online) The transmission vs. signal amplitude for the loss- 
compensated metamaterial of a three-layer system (NIM - gain material - NIM) with 
gain bandwidth of 5THz, for different pumping rates r pump . r pump is increased from 
4.0 x 10 9 s" 1 (lowest) to 5.5 x 10 9 s~ 1 (highest) in steps of 1.0 x 10 8 s _1 . The material 
parameters for NIM are s = /i = — 1 + 2i. The metamaterial response is linear in a 
very wide range. When the loss-compensated transmission reaches exactly unity, the 
pumping rate r pump = 4.70 x 10 9 s~ 1 , which is called the critical pumping rate. For 
incident fields stronger than 10 V/m the metamaterial behaves nonlinearly. 

signal is only affected by the dispersion but not dependent on the signal strength. For 
a three- layer system (NIM - gain material - NIM), we have calculated the transmission 
versus the strength of the electric field of the incident signal for several pumping rates 
close to the critical pumping rate r pump = 4.70 x 10 9 s _1 , as shown in figure [TOl We 
found it has a linear response within a very broad range up to incident electric field of 
10 3 V/m. If we use 19 layers of figureEl the critical pumping rate is 1.98 x 10 9 s _1 , which 
is even smaller than half of the three-layer case, and the linear regime becomes narrower 
and drops faster than the 3-layer case for higher strength of incident electric field (shown 
in figure [TTT) . To include the nonlinearity of gain material for strong incident signal, it 
is necessary to do a self-consistent calculation using FDTD method. 

As an example, we have also studied the three-layer system with different losses, to 
see how much r pump we need to compensate the losses. Figure [12] shows there exists a 
linear relation between the critical pumping rate and the imaginary part of the refractive 
index n of NIMs. 

We have also numerically calculated the susceptibilities of the gain material to 
see if it really has a Lorentzian lineshape. We first let a Gaussian pulse of a given 
amplitude (10 V/m) propagate through the metamaterial system and calculate the 
time-domain electric polarization P(r, t) and the local electric field E(r, t). Then 
we implement the Fourier transforms to obtain the frequency-domain polarization 
and electric fields and calculate the frequency-dependent susceptibility by using the 
equation ~x!{lS) + ix"^) — Pi^) / £ qE(uj). Simulations are done for both 3 and 19 
layers and results are compared with the analytic results calculated using the equations 
p3] x > = -x»Ax/(l + Ax 2 ) and x" = Xo/i 1 + A % 2 ) with Ax = 2 - u a )/T a and 
Xq = —cr a AN/(e u> a r a ), where AiV = N 2 — N%. As shown in figure [T3| we found that 
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Figure 11. (Color online) The transmission vs. signal amplitude for the loss- 
compensated metamaterial of figure [5] with gain bandwidth of 5 THz at the critical 
pumping rates r pump = 1.98 x 10 9 s -1 . The material parameters for NIM are same 
as figure ITD1 For incident fields stronger than 10 3 V/m this metamaterial becomes 
non-linear. 
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Figure 12. (Color online) The critical pumping rates for different imaginary parts 
of the refractive index n of NIMs. The structure is a three-layer system (NIM - gain 
material - NIM) and s = fi for NIMs. 



the numerical susceptibilities are the same as the analytic ones and they do have a 
Lorentzian lineshape. 

To understand the lasing behavior of the metamaterial system, we increase the 
pumping rate to provide more gain from the gain material. We found the amplification 
of the incident signal gets larger and at last the system starts lasing till the pumping 
rate reaches a certain high value. Figure 14(a) shows the lasing behavior at the pumping 
rate r pump = 1.5 x 10 10 s _1 and there appears a peak at the emission frequency 100 THz 
in the corresponding Fourier transform (figure 14(b)). 

If we use real metal layers instead of negative index materials, we can not 
compensate the losses of the metals. The reason is that the permittivity e for metals is 
large and negative and we'll have large reflections due to the impedance mismatch. 
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Figure 13. (Color online) The numerical and analytical results for the susceptibilities 
of gain materials as a function of frequency, (a) a three-layer system (NIM - gain 
material -NIM) at the critical pumping rate r pump = 4.7 x 10 9 s _1 . (b) a 19-layer 
system of figure [9] at the critical pumping rate r pump = 1.98 x 10 9 s _1 . 




(a) (b) 

Figure 14. (Color online) (a) the time-dependent electric field of the transmitted 
wave for a three-layer system (NIM - gain material - NIM) and (b) the corresponding 
Fourier transform in frequency domain for the lasing in (a). The pumping rate 
Tpump — 1-5 x 10 s . 



3.3. One layer of gain material embedded in a square lattice of Lorentz dielectric 
cylinders 

In section 13721 we simply force the permittivity and the permeability of the metamaterial 
to be negative to have an unrealistic negative index material. In this section, we consider 
a two-dimensional (2D) metamaterial system (shown in figure fT5l) which consists of one 
layer of gain material and two layers of dielectric wires that have a Lorentz-type resonant 
electric response to emulate the resonant elements in a realistic metamaterial, such as 
cut- wires. 

We first study the three-layer system of figure [15] with the full width at half 
maximum (FWHM) 5THz for Lorentz dielectric (i.e., / = 2.5 THz) and 20THz for 
gain. The transmission, T, the reflection, R, and the absorption, A = 1 — T — R, as a 
function of frequency for the system are obtained in the propagation direction. With the 
introduction of gain, the absorption at the resonance frequency of 100 THz decreases 
and reaches at a certain pumping rate. So the gain compensates the losses. If we 
continue increasing the gain, i.e., the pumping rate, the system gets overcompensated 
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Figure 15. (Color online) One layer of gain material (orange) embedded in a square 
lattice of dielectric square cylinders (blue) that have a Lorentz behavior. The dielectric 
constant of the cylinders is given by e = 1 + u>p/(u>p — 2iujj — uj 2 ), where the resonance 
frequency f p — lu p /2tt — 100 THz and 7 = 2tt/, and / takes different values in the 
cases we have examined. The dimensions are a = 80 ran, Wl = 40 nm, and w g = 30 run. 
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Figure 16. (Color online) (a) The retrieved results for the real and the imaginary 
parts of the effective permittivity e with gain and without gain. Below compensation, 
t = 0.89; gain and Lorentz bandwidths are 20 THz and 5 THz, respectively, (b) and (c) 
are same as (a), but for the loss compensated case (t = 1) and the overcompensated 
case (t = 1.34), respectively, (d) - (f) are same as (a) - (c), respectively, except with 
gain bandwidth 5 THz and Lorentz bandwidth 20 THz. 



and the absorption becomes negative. To see how the losses of the emulated resonators 
get compensated by the gain, we exploit the usual retrieval procedure based on inverting 
the scattering amplitudes p£] to obtain the effective permittivities e without gain and 
with gain. Figure [ToT a) - 16(c) show the retrieved results for the real and the imaginary 
parts of the effective permittivities e of the system for the below compensation, loss- 
compensated, and overcompensated cases, respectively, together with the effective 
permittivity without gain. The retrieved results for e without gain have exactly the 
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Lorentzian shape but the amplitude of the real and the imaginary parts of e is a factor 
of 4 less than the Lorentz formula for the square cylinders. This is due to the filling 
ratio of the square cylinders in the unit cell. Due to the loss compensation from the gain 
material, one can see the imaginary part of the effective permittivity gets lower as the 
gain increases. Below compensation, its value at the resonance frequency of 100 THz 
is positive, while it's zero and negative for the loss-compensated and overcompensated 
cases, respectively. Notice that we can have Re(e) ~ 2 with Im(e) ~ at 98 THz and 
Re(e) ~ —1 with Im(e) « at 101 THz for the below compensation case (figure 16(a)). 
For the loss-compensated case (figure 16(b)), we have Re(e) ~ 1 with Im(£) ~ at 
the resonance frequency 100 THz, which is just same as vacuum and makes no sense 
for us. For the overcompensated case (figure 16(c)), the imaginary part of the effective 
permittivity e is negative within all the frequency range. 

Second, we study the three-layer system with the FWHM of gain smaller than 
Lorentz dielectric, where the bandwidths for gain and Lorentz dielectric are 5 THz and 
20 THz (i.e., / = 10 THz), respectively. The introduction of gain develops a peak at 
the resonance frequency of 100 THz for the transmission while the absorption has a dip. 
The retrieved results for the real and the imaginary parts of the effective permittivities 
e without gain and with gain are plotted in figure 16(d) - 16(f) for the three different 
cases discussed above. Similar to the first case we examined where the loss bandwidth is 
smaller than the gain, the imaginary part of the effective permittivity e gets smaller due 
to the gain. The difference is we get interesting results for the overcompensated case 
instead of the below compensation case, where we can have Re(e) ~ with Im(e) ~ 
at 97 THz and Re(e) « 2.1 with Im(e) « at 103 THz. 

So for both the two systems, one can obtain a lossless metamaterial with positive 
or negative Re(e), either below compensation or over compensation. In figure[16]we also 
have plotted the sum of Im(e) without gain and the imaginary part of e g , the dielectric 
function of the gain material. One can see the imaginary part of e of our total system 
with gain is equal to the sum of Im(e) and lm(e g ). This is unexpected because there 
is no coupling between the Lorentz dielectric and the gain. This is indeed true for the 
2D Lorentz dielectric cylinders, because they have a continuous shape like a solenoid 
and the gain material slabs have zero depolarization field. Different from finite length 
wires [hence a three-dimensional (3D) problem] where the dipole interactions between 
Lorentz wires and gain material are dominated by the near field 0(l/r 3 ), the interaction 
for infinite length wires is only via the propagating field 0(u\n \kr\), and much weaker. 
That's why the Lorentz wires and the gain material are approximately independent in 
our 2D simulations. So there is a need for a true 3D simulation to solve this problem and 
obtain different behaviors. However, the 3D simulation is computationally excessively 
demanding. 

Like the layered system in section 13.21 if we keep increasing the pumping rate, 
i.e., the gain, at last both of the two systems will have lasing. For example, when the 
pumping rate reaches r pump = 3.2 x 10 10 s _1 , the three-layer system with gain bandwidth 
5 THz and Lorentz bandwidth 20 THz starts lasing — the system itself has a coherent 
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Figure 17. (Color online) Geometry for a unit cell of the square SRR system with gain 
embedded in the gap (orange). The dimensions are a = 100 nm, I = 80 nm, t = 5nm, 
d = 5 nm and w = 15 nm. 



self-sustained steady output. 

3.4- 2D split ring resonators (SRRs) with gain material inclusions 

To avoid the decoupling problem in section 13.31 and still limit our simulations in a 2D 
model, we consider a 2D split ring resonator (SRR) as a more realistic and also more 
relevant model, where the relevant polarization is across the finite SRR gap and therefore 
the coupling to the gain material is dipolelike. In figure ITTJ, we show the geometry for 
the unit cell of the SRR system with the gain material embedded in the gap. The 
dimensions of the SRR are chosen such that the magnetic resonance frequency of the 
SRR overlaps with the emission frequency (100 THz) of the gain material. Due to the 
strong electric field inside the gap, there will be strong coupling between the SRR and 
the gain material. We also want to see if the losses of the magnetic response can be 
compensated by the electric gain. 

Simulations are done for one layer of the square SRR with gain bandwidth of 



20 THz. Figure 18(a) shows the retrieved results for the real and the imaginary parts 
of the effective permeability fi, with gain (pumping rate r pump = 1.0 x 10 9 s~ 1 ) and 
without gain. One can see that with the introduction of gain, the gain undamps the 
magnetic resonance of the SRR and the weak and broad resonant effective permeability 
\i of the lossy SRR becomes strong and narrow. The FWHM with gain is 2.61 THz, 
while the FWHM without gain is 5.85 THz, which is more than twice larger than the 



former. Notice that in the off-resonance range in figure 18(a) , we can obtain the effective 



permeability fi with a smaller imaginary part with the introduction of the gain, which 



means the magnetic loss is compensated by the electric gain. Figure 18(b) shows the 
retrieved results for the real and the imaginary parts of the corresponding effective 
index of refraction n, with and without gain. Note that for a lossless SRR, n is purely 
real away from the resonance except in a small band above the resonance where it's 
purely imaginary due to the negative \i. At the frequency of 96 THz, slightly below the 



resonance (figure 18(b)), the imaginary parts of the index of refraction n without and 
with gain are 1.36 and 0.754. Then we can find the effective extinction coefficient without 
gain is a = (u/c)lm(n) pa 2.74 x 10 4 cm _1 and the one with gain is a pa 1.52 x 10 cm" 1 . 
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Figure 18. (Color online) The retrieved results for the real and the imaginary parts of 
(a) the effective permeability \x and (b) the corresponding effective index of refraction 
n, with and without gain for a pumping rate r pump = 1.0 x 10 9 s" 1 and the SRR system 
of figure |T7l The gain bandwidth is 20 THz. Notice that the width of the magnetic 
resonance with gain is 2.61 THz. 



And hence the effective amplification coefficient of the gain in the combined system is 
a ~ —1.22 x 10 4 cm -1 , which is much larger than the amplification a ~ —9.2 x 10 2 cm" 1 
for the bulk gain material [H] at the given pumping rate r pump = 1.0 x 10 9 s _1 . This 
is due to the strong local electric field enhancement in the gap of the resonant SRR. 
While we have the incident electric field 10 V/m, the induced electric field in the gap is 
around 450 V/m. Taking the observed field enhancement factor around 45 in the gap of 
SRR, the energy produced by the gain in the gap is around 12 times larger than by the 
homogeneous bulk gain material in the size of the unit cell, which agrees well with the 
factor around 15 between the simulated SRR effective medium and the homogeneous 
gain medium. If we continue increasing the pumping rate, the magnetic resonance 
becomes narrower (0.96 THz for pumping rateT pump = 1.8 xlO 9 s" 1 ). When the pumping 
rate reaches r pump = 1.9 x 10 9 s _1 , the metamaterial system gets overcompensated 
and the imaginary part of the effective permeability /i at the resonance frequency gets 
flipped down and becomes negative. If we increase the pumping rate even more (around 
r pU mp = 5.0 x 10 9 s" 1 ), the SRR system starts lasing [191120] . 

Instead of having the gain in the gap of SRR, we also have done simulations on the 
SRR system where the gain is surrounding around the SRR. The retrieved results for 
the real and the imaginary parts of the effective permeability // without gain and with 
gain for a given pumping rate T pump = 8.0 x 10 8 s _1 are plotted in figure [191 Due to 
the larger gain filling in the unit cell, the coupling between the gain material and the 
SRR gets stronger and the losses of SRR are easier to be compensated. Compared with 
the case with the gain in the gap only, the effective permeability \x of the simulated 
SRR system in figure [T9l can have stronger and narrower resonance (FWHM = 1.6 THz) 
even with lower pumping rate. Similar to the previous case, the magnetic resonance 
becomes stronger and narrower (0.7 THz at r pump = 1.0 x 10 9 s _1 ) as we increase the 
pumping rate, and the system gets overcompensated so that the magnetic resonance 
peak gets flipped down when the pumping rate reaches r pump = 1.1 x 10 9 s _1 . It's also 
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Figure 19. (Color online) The retrieved results for the real and the imaginary 
parts of the effective permeability fi with and without gain for a pumping rate 
Tpump = 8.0 x 10 8 s -1 and the SRR system where the SRR is surrounded by gain. 
The gain bandwidth is 20 THz. Notice that the width of the magnetic resonance with 



gain is 1.60 THz, narrower than figure 18(a) 



easier for this system to have lasing [19j [20] , which is observed when the pumping rate 
is ~ Tpunip —— 2.5 x X0 s 

4. Conclusions 

We have performed numerical simulations on one dimensional gain material system 
by using FDTD method. The system starts lasing when the input reaches over the 
lasing threshold. For more input power and wider gain slab, the lasing is faster. 
Comparisons were done for a gain material slab between the optical pump method and 
its homogeneous pumping rate simplification and results show that this simplification 
can be valid for a thin gain slab. 

A self-consistent model incorporating the gain into dispersive metamaterial 
nanostructure was proposed and numerically solved. We numerically show that the 
losses of dispersive metamaterials can be compensated by gain by investigating the 
transmission, reflection and absorption data as well as the retrieved effective parameters. 
There is a relatively wide range of input signal amplitudes where the metamaterial-gain 
system behaves linearly. When the amplitudes get higher, the system becomes nonlinear 
due to the nonlinearity of the gain material itself. It's necessary to have self-consistent 
calculations to determine the signal range where we can expect linear response. Further, 
if we have strong signals so that we are in nonlinear regime or we want to study lasing, 
self-consistent calculation is needed. As examples, two SRR systems with different gain 
inclusions were studied. We have demonstrated that the magnetic losses of the SRRs can 
be easily compensated by the electric gain. The pumping rate needed to compensate the 
losses is much smaller than the bulk gain material. The losses of the SRR surrounded 
by gain can be easier to be compensated than the SRR with gain in the gap only due 
to more coupling with the gain. Provided that the pumping rate is high enough, the 
metamaterial nanostructures can have lasing. 
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